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Abstract 

We discuss the thermal entanglement close to a quantum phase transition by analyzing the 
pairwise entanglement for one dimensional models in the quantum Ising universality class. We 
demonstrate that the entanglement sensitivity to thermal and to quantum fluctuations obeys uni- 
versal T 7^ 0-scaling laws. We show that the entanglement exhibits a peculiar universal crossover 
behaviour , shedding light on the mechanisms bringing quantum criticality up to finite tempera- 
tures. 
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Quantum phase transitions (QPTs) are zero temperature transitions driven by quan- 
tum fluctuations Even though absolute zero is not accessible experimentally, QPTs Jiave 



profound influences at flnite temperature, making the phenomenon detectable 
It is the relative role of the thermal over quantum fluctuations inherited from T = 
that constitutes a characterization of the system at low temperature^]. The interplay be- 
tween classical statistical mechanics and quantum mechanics may be so tight that a qual- 
itatively new phenomenology can emerge. Thereby it is expected that the entanglement 
could play an important role, comp 
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ementing traditional approaches employed in statisti- 
Q, ll^ • Though entanglement is generically present 



in quantum many particle systems, the key to distill it quantitatively out of the correla- 
tions was provided only recently setting the so called "entanglement measures" , assigning a 
precise number to the entanglement encoded in a given state 
glement in systems of many quantum particles now constitutes a new angle of research in 



16 



171]. Analyzing the entan- 



the physical literature 18|, Il9|]. In the extreme case of T = it was shown that entanglement 
and correlation functions share the property to be sensitive to the qualitative change of the 
ground state at the quantum critical point, and with the same non-analytic behaviour; how- 
ever, two-points entanglement and correlation functions are distinct in that the flrst can be 
short-ranged whenever the system is correlated to all scales}?]. Many aspects of QPT were 
investigated under this light. Various QPT in spin- 1/2 models have been studied 



10, 



111; 



the string order parameter in gaped spin 1 systems has been interpreted in terms of lo- 
calizable entanglement 12^. Connection of the entanglement entropy with Conformal Field 

nn 

Theory has been formulated [1 31. 1 141] . Finally entanglement in QPT of itinerant electrons was 
proved to be related to ground state charge and spin susceptibilities[l5|]. 

Studies on entanglement close to QPT are conflned to zero temperature. The aim of 
the present study is to explore quantitatively the role of the temperature on the quantum 
critical properties of entanglement encoded in the state of the system. We refer to the flnite 
temperature entanglement and not to ordinary correlation functions (containing both clas- 
sical and quantum correlations) to describe the effects fanning out from the quantum phase 
transition, that is a pure quantum phenomenon. We demonstrate that the entanglement 
sensitivity to thermal and to quantum fluctuations obeys universal T ^ 0-scaling laws. We 
show that the entanglement, together with its criticality, exhibits a crossover behaviour with 
new aspects of it, which have universal features within the class of models we considered. 
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Finally we elaborate on how typical is the crossover for the entanglement. 



Concretely, we focus on a spin chain whose Hamiltonian is[20|] 



N 



N 



H = - J ^(1 + 7)5r^^i + (1 - l)S!Sl, - h J2 s^- 



(1) 



1=1 



i=l 



The spins 5"° 



0" , a 



X, y, z {a" are Pauli matrices) experience an exchange interaction 



with coupling J > and uniform magnetic field of strength h. Besides to model actual 
compounds 211] , Eq.([l]) constitutes a standard arena in statistical mechanics for reliable 
studies in the theory of quantum critical phenomena. In fact at T = 0, for N ^ oo the 
systems ([1]) have a quantum critical point at a = h/2J = ac = 1, and for < 7 < 1 they 
identify the quantum-Ising universality class. The low temperature behavior is influenced by 
the interplay between thermal and quantum fluctuations of the order parameter (S^). The 
renormalized-classical crosses-over the quantum disordered regimes through the so called 
"quantum critical region" [l|, ^ In the T — a plane a l^-shaped phase diagram emerges, 
characterized by the cross-over temperature customarily defined as T^ross = la — ad^'^ (the 
critical exponents z/ = 2; = 1 for the model ([1])) fixing the energy scale [23]. For T ^ 
Tcross the thermal De Broglie length Xth is much smaller than the average spacing of the 
excitations C,c', therefore the correlation functions factorize in two contributions coming from 
quantum and thermal fluctuations separately. The quantum critical region is characterized 
by T ^ Tcross- In there Xth ~ ^c, and the factorization of the correlation functions does not 
occur. In this regime the interplay between quantum and thermal effects is the dominant 
phenomenon affecting the physical behaviour of the system. Here we describe how such 
interplay is manifested by the finite temperature two-spins entanglement. We consider the 
entanglement between two spins, at chain positions i and j. This can be expressed by 
the concurrence C{R) = max{0, Ao(-R) — Ai(i?) — A2(-R) — A3(-R)}, where A„ are the 
square roots of the eigenvalues of the matrix ay ® ayp*ay (g) ayp, in descending order, being 
p the reduced density matrix and R = \i — j\. We point out that for the models ([1]), 
C{R) = C{R, a, T) can be obtained analytically at any T and a since p is accessible exactly 
by using results of the Refs.[26|, l27|. For generic 7's, close enough to the quantum critical 
point, the range of thermal concurrence results to be constant; deviations are ultimately 
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a/1 — 7^ that turns out a non 
29I . We remark that the customary analysis |l, [g, [22I based on the 



due to the factorization of the ground state at a 
critical effect lol . 



correlation length and the asymptotics of the correlation functions is not feasible here, since 



the entanglement extends over a short range, that is found non-universal{3]. Despite of this, 
we shall evidence how the crossover is reflected at short-distance, on the entanglement. 

The first aspect we consider is the scaling of the entanglement close to the quantum 
critical point. It was demonstrated that the trait of the QPT can be captured by the field 
derivatives daC{R) diverging logarithmically at T=0 and a a^, with scaling behaviour [3]. 
At finite temperature we study both daC and OtC. Following the standard paradigm (see f.i. 
Refs.jl.l30l) the parameters a and T are indicative for the 'strenght' of thermal and quantum 
fluctuations of the order parameter respectively. Then daC and dxC can be interpreted as 
the response of the entanglement when those fluctuations are tuned. Here we analyze the 
scaling properties of such two quantities in a finite temperature neighborhood of the quantum 
critical point. We first consider the Ising model (7 = 1). Also at T 7^ C(l) and C(2) are 
the non vanishing concurrence. 

For daC{R), at T 7^ the divergence at the critical point is rounded off to maxima 
scaling as ~ Oc — T^'^^, that are the finite temperature effects remaining out of the QPT. 
According to the theory of the quantum critical phenomena j^, Q, 24 1, the scaling ansatz is 

daC{l) ^ ln[T^g(-^)] (2) 

cross 

where T is a non universal exponent; Eq.(l2]) was verified in Fig{T] where the data collapse of 
Q as function of T /T^ross is presented. Universality in the critical entanglement was verified 
by analyzing ([2]) for different 7's (the results are not shown here). C(2) shows the same 
universal behaviour. 

Though dxC^R) is not diverging at the QPT it is affected by the quantum criticality. In 
fact we observe that 

drCil) ^ P(-^) (3) 

cross 

The scaling function \P{T/Tcross)\ is characterized by pronounced maxima at certain 
crossover temperature T* = aTcross (see the caption of Fig 12]). Close enough to the quantum 
critical point the position of the maxima is found to be independent of 7 and R. There- 
fore, by gaiting into the quantum critical region the entanglement crosses-over at T = T*, 
experiencing the largest variation for thermal fluctuations. This effectively "separates" the 
entanglement in the renormalized classical and disordered regimes from that one in the 
quantum critical regions. However, as we shall point out below, such separation is much 
more involved. 
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To study the symmetry of the entanglement in the a — T plane we analyzed the directional 
derivative D^C^R) = |VC(i?) ■ u| = IdrC^R) sin a + daC{R) cosa| as function of the slopes 
u = (cos a, sin a), tana = T/Tcross emerging from the quantum critical point (see the 
caption of FiglS]). The behavior of daC{R) and dTC{R) implies that D^^C{R) is very small 
along the critical line a = ac, within a cusp-like domain. In a sense within such a region, 
close to the quantum critical point, the pairwise entanglement is particularly "rigid". This 
could provide a mechanism explaining how quantum correlations are robust up to finite 
temperatures particularly along slopes within the quantum critical region. 

Now we analyze how the entanglement is affected by the interplay between quantum 
and thermal effects. We consider the quantity: drC (R) / daC (R) (we note that a formally 
similar quantitity, the so called Grueneisen parameter, emerged in recent studies of quantum 
criticality 25|). This quantity can be interpreted as a measure of the relative role of the 



effects of thermal and quantum fluctuations of the order parameter on the entanglement. 
The analysis evidences two regimes (FigHj). daC{R) is dominant at very low temperature 
and near the critical line a = ac- Instead in the crossover to the quantum disordered and 
renormalized classical regimes the role of entanglement's thermal sensitivity is enhanced. 
The importance of the effect of temperature in the entanglement crossover can be pointed 
out also analyzing how the critical divergence of daC{R) is affected by variations of the 
temperature: dT[daC{Ry\. Due to the vanishing of the gap at the quantum critical point, 
in the region T ^ Tcross an arbitrarily small temperature is immediately effective in the 
system. Interestingly enough we observe that dT[daC{R)]is enhanced also along the line 
Tm = PTcross with (3 < a and spread out on a finite range of T /Tcross- Td < T < Tc2 (see 
the caption of FigJSj); the values of Tm, Td and Tc2 are proved to be independent of 7 and R- 
We notice that Tm characterizes also the crossover in FigJUthat is approximately enclosed 
between T* and Tm- The analysis of the signs of dT[daC{R)] leads to the conclusion that 
for T > T* a reduction of the temperature corresponds to an increase of the effects of the 
critical behavior of daC{R)] while for T < T* a reduction of the temperature effectively leads 
the system away from the criticality, consistently with the important role of the energy gap 
scale in such a regime |l, For T < Td the temperature does not affect daC{l). 

To estimate the impact of the entanglement in the total correlation we apply a variation 
of the arguments developed in {32]. Specifically, we analyze the quantum mutual information 
amoung the i — th and j — th spin: lij = Si + Sj — Sij, with Sx = Trp^ logpx being the 
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Von Neumann entropy of px 3l| , as a valuable measure of quantum & classical correlations, 
taking into account of all two-point spin correlation functions. It emerges that a similar 
phenomenology found for the concurrence is experienced by I if, moreover we found that 
the crossover temperatures for lij are shifted by less than 10% respect to the corresponding 
temperatures for the the entanglement. Then: By distilling out the quantum part from the 
total correlation, it emerges how the crossover is little affected by the classical correlations, 
suggesting that it is driven solely by the interplay between thermal and quantum sensitivities 
of the entanglement . 

The present analysis reveals that the entanglement sensitivities to thermal and quantum 
fluctuations close to a QPT are universal functions of the angle at which the critical point is 
approached. By studying how the entanglement in the ground state of the system is affected 
by the contribution coming from the excitations, we analyzed how pure quantum effects are 
affected by the finite temperature. In particular, it is shown that the entanglement exhibits 
a crossover behaviour evidencing a "fine structure", characterized by energies T* and Tm 
that are found independent of the specific model in the class (IH). Interestingly enough, a 
similar anomaly in the entanglement sensitivity was evidenced recently in the context of the 
spin-boson model. In fact it was proved that the coherent-incoherent oscillations crossover 
is accompained by the largest variation of the entanglement between the two-level system 



and the bosonic bath 



33| (see also [3J]) 



Our scaling result together with the study of the entanglement crossover are particularly 
significant, in our opinion, in that they involve pure quantum correlations at a short-range; 
this is pointed out by analyzing also the quantum mutual information. /^From our analysis a 
neat meaning is sheded on the specific crossover temperature we found, beyond its standard 
interpretation as a mere reference energy scale for the crossover. 

In statistical mechanics, our analysis can constitute a crucial ingredient for a deeper 
understanding of the role of the quantum distilled part of the correlations entanglement 
in the crossover phenomena close to QPTs. As drawn in jsS, 36], spin-off in the research 
on the mechanisms behind the puzzling low temperature behavior of exotic materials like 
high Tc superconductors and heavy-fermions are possible. We point out that, being at 
finite temperature, our study can open the way towards an experimental observation of the 
physical effects that entanglement has on the criticality of the system. We believe that 
the present analysis could be important both for pure research and quantum computation 
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perspectives. 
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FIG. 1: Scaling of entanglement at finite temperature close to the quantum critical point. The 
scaling function Q{T /Tcross) for daC{l) is presented. Temperatures are expressed in units of J. 
Q = exp{daC {!)) /T'^ as function of the non-scaled coordinate T is shown in the insets; different 



colors correspond to different distances d = (a — Oc)^ + from the quantum critical point. 
Plotted data belong to a neighborhood of the QCP of radius 10^^ < d < 10^^ (respectively violet- 
black data); data closer to the critical point collapse into the same line asymptotically, thus proving 
the scaling ansatz (the deviations are corrections to the scaling due to inspections at larger distances 
to the critical point). The functional form of Q is obtained by imposing that daC{l) reduces to 
the T = resultQl; for that: Q{T 0) ^ a{T /Tcross)''^ ] Eq.© reduces to Aln\a - ad + B, 
fixing B = — Ina, A = zi/T; v e z are related to T at a = Oc since da\a=acC{l) = TlnT -|- cost, 
implying that zi^ = A/T. We found da\a=a,C{l) = -0.2701 ln(r) - 0.2089 consistently with the 
critical exponents and the results of [3] ■ 
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0.1 1 10 10 1 0.1 



T/T T/T 

cross cross 

FIG. 2: The scaling properties of dxC^l) (7 = 1) are verified for 10~^ < < 10^^ around the 
quantum critical point. Such properties are found to be universal and independent of R. The 
middle inset evidences that the maxima are independent of 7 ; the corresponding temperatures are 
plotted as functions of the distance to the critical point d, for different values of 7. At small d the 
maximum position saturates at T* = aTcross with a ~ 0.595 ± 0.005. 
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FIG. 3: Symmetry of the entanglement in the T—a plane. The directional derivative DuC{R), with 
tana = T/Tcross- We observe that the critical point can be approached with constant DuC{R) 
following cusp-like trajectories. The density plot refers to the Ising model 7 = 1 and to nearest 
neighbour concurrence C(l). T = T* and T = Tm are drawn as dashed and thick lines respectively. 
The same behaviour of the entanglement arises at any values of 7 and R. 
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FIG. 4: The relative role of thermal over quantum fluctuations of the order parameter on the 
(9tC(1) 

entanglement for 7 = 1. The crossover regions between the low temperature and quantum 

C'aL' (1) 

critical regions are approximately enclosed in Tm <T* which are drawn as thick and dashed 
lines respectively. We notice that such a ratio has also a geometrical meaning in that it can be 
obtained as the angle of steepest concurrence variation: da[D^C{R)] = 0. Therefore the two 
regions where daC{l) is dominant correspond to abrupt concurrence changes at = 0. line). 
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FIG. 5: The effect of temperature on the anomalies originated from the critical divergence of the 
field-derivative of C{R) can be measured by dT[daC{R)]. The density plot corresponds to 7 = 1 
and i? = 1. T = T* and T = Tm are drawn as dashed and thick lines respectively. Maxima below 
T* are found at Tm = f3Tcross with (3 ~ 0.290 it 0.005 and they are independent of 7 and R; the 
crossover behaviour is enclosed in between the two flexes of dT[daC{R)] at Td Tc2'-, such values are 
fixed to: Td = (0.170 ± 0.005)Tcross and Tc2 = (0.442 it 0.005)Tcross and found to be independent 
of 7 and R. For T < T^i dT[daC{R)\ ~ 0. Scaling properties are inherited in dT[daC{R)\ from 
daC{R)] (see Eq.©). 
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